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<D ■ ABSTRACT 

For students familiar with equilibrium statistical mechanics, the notion of a positive specific heat, 
being intimately related to the idea of stability, is both intuitively reasonable and mathematically 
provable. However, for system in non- equilibrium stationary states, coupled to more than one energy 
C/3 ' reservoir (e.g., thermal bath), negative specific heat is entirely possible. In this paper, we present a 

"minimal" system displaying this phenomenon. Being in contact with two thermal baths at different 
temperatures, the (internal) energy of this system may increase when a thermostat is turned down. 
In another context, a similar phenomenon is negative conductivity, where a current may increase 
by decreasing the drive (e.g., an external electric field). The counter-intuitive behavior in both 
processes may be described as " getting more from pushing less." The crucial ingredients for this 
phenomenon and the elements needed for a "minimal" system are also presented. 
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I. INTRODUCTION 



> 

■ To most physicists, it is intuitively reasonable that the average internal energy, U, of a system in contact with 
a thermal bath should increase with T, the bath temperature. For those well versed with equilibrium statistical 
I/"") , mechanics, it is trivial to prove this statement, usually in the form of the specific heat, dll/dT, being positive. 
Indeed, this notion is so deeply rooted in the idea of stability that its contrary may be rejected immediately, even for 
non- equilibrium steady states. In this paper, we show the presence of a negative specific heat in a minimal system, 
coupled to two thermal baths at different temperatures. We also point out a crucial ingredient which makes this 
phenomenon possible, namely, the presence of "barriers" which trap the system in a higher energy state. In a similar 
vein, we discuss negative conductivity, using again a simple model to demonstrate the phenomenon of increasing 

S current by decreasing the drive (e.g., an external electric field). Here too, a crucial ingredient is the presence of 
"obstructions" which trap the system in a low current state. Noting that both processes involve a similar idea, we 
'-^J ■ coin a "folksy" phrase - getting more from pushing less - to describe this class of counter-intuitive behavior. 

We begin with a brief summary of standard equilibrium statistical mechanics, for which no dynamics is needed. 
In contrast, for non- equilibrium steady states, the underlying dynamics (which leads a system into such states) is 
crucial. The most convenient framework for describing both scenarios is the master equation , which we review in the 
next section, for the convenience of those readers unfamiliar with this approach. In this context, a clear distinction 
between equilibrium and non-equilibrium steady states will be presented. In section 3, we display the presence of 
negative specific heat in some minimal systems and a more typical, many-body system, which are coupled to two 



H 

thermal baths. Similarly, negative conductivity is shown to arise in a simple hopping model in Section 4. Concluding 
in the last section, we recapitulate the key ingredients needed for such phenomena and draw an intuitive picture of the 
mechanism. For ease of reading, the details of the algebra are taken out of the main text and placed in an appendix. 

II. CONTRASTS BETWEEN EQUILIBRIUM AND NON-EQUILIBRIUM STEADY STATES 

To describe a physical system statistically, a complete and relevant set of microstates (or configurations: C) for the 
system must first be specified. For simplicity, let us consider finite systems with finite degrees of freedom (e.g., finite 
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Ising models), so that the C's form a discrete, countable set. With appropriate mathematical tools, our considerations 
can be extended to the thermodynamic limit of systems with continuous degrees of freedom. Next, we need P[C], the 
probability for finding the system in each C. From here, the average of an observable quantity (e.g., internal energy), 
which takes the value 0[C] in configuration C, is then given by 

{0) = J20[C]P[C}. (1) 

c 

Beyond averages, fluctuations and correlations are similarly computed. 

When a system is evolving, then the probability distribution is a function of time: P[C; t], leading to time dependent 
averages, fluctuations, etc. Clearly, an "equation of motion" is needed to describe the evolution of P[C; t]. In principle, 
we should derive such an equation from, e.g., Newton's equations. In practice, however, this is typically not feasible. 
Instead, let us follow the standard framework of a master equation [jD: 

|P[C; (C - C) P[C; i] - JJ (C - C) P[C; *]} (2) 

which is basically a continuity equation for the probability density. Here, R (C — > C) stands for the rate a configuration 
C changes to C' . Most frequently, these rates are not time dependent, i.e., the dynamics is invariant under time- 
translation. As in the case of Eq. (|^) itself, they can be found in principle, once we specify how our system is coupled 
to its environment (e.g., a thermal reservoir). Again, in practice, the task of finding these rates is prohibitively 
complex, so that progress in this approach relies, typically, on postulating reasonable forms, based on sound physical 
principles. Returning to Eq. (|^), we see that, being linear, the right hand side may be written as a matrix, L, 
operating on a "vector." In short hand (with \P) t = P[C;t]), this equation takes the form 

d t \P) t = L\P) t . (3) 

Explicitly, the matrix elements arc 

Lc,c =R(C'^C)-6 (C, C) Y, R ( C ^ c ") ( 4 ) 

C" 

where 6 (C,C) is the Kronecker delta. 

For simplicity, let us focus on systems which eventually settle down into a unique time- independent state, i.e., 



Urn P(C,t) =P*{C). (5) 



In other words, P* satisfies dP* /dt = and is the eigenvector of L with zero eigenvalue. Further, the real parts of 
all other eigenvalues are strictly negative, so that P* is a stable state. 

Clearly, this approach should embrace systems which evolve towards thermal equilibrium. If we wish to reproduce 
the results from the theory of equilibrium statistical mechanics (e.g., the Boltzmann factor in the canonical ensemble), 
there must be some constraints on these rates. These constraints are known as detailed balance. Specifically, denoting 
the energy of a configuration of our system be given by the Hamiltonian W[C], we must impose, for an isolated system 
with total energy E, 

R(C -» C) = R{C -»■ C) (6) 

if H[C') = H[C] = E and R(C — > C) =0 otherwise. Then P*(C) oc 1, which is the fundamental hypothesis (micro- 
canonical ensemble), clearly satisfies L|P*) = 0. Generalizing to the canonical case (with (3 = l/Zc^T), we demand 

^^l = e X p{p(H[C'}-H[C})} (7) 

and verify that a distribution given by the Boltzmann factor 

P: q (C) oce-™ (8) 

indeed satisfies L \ Peq) = 0- Apart from these constraints, there is much leeway in postulating the i?'s, corresponding 
to the fact that equilibrium states are independent of the details of the dynamics. In particular, it is not crucial that 
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R (C — ► C) 7^ for all pairs of (C, C), though there must be enough non-zero rates so that every configuration (within 
the desired ensemble) may be reached from any other one. As an explicit example, we give the Metropolis rate 0: 

J2(C — C) = min[l, ^(WPI-WM)] , (9) 

which is the basis of numerous successful Monte Carlo simulations of systems in equilibrium. 

From the condition of detailed balance, a stronger statement about an equilibrium state emerges. Not only is the 
distribution time- independent, the net (probability) current between any pair (C,C) vanishes! To be more explicit, 
we regard the right hand side of the master equation (||) as a sum of the net currents (from C into C ) : 

R (C -> C) P[C; t]-R(C^ C) P[C; t]. (10) 

Then detailed balance, embodied in Eqs. (|^) or (0), implies the vanishing of all net currents. An analog of this 
situation in electrodynamics is electrosfaiics, where the charge distribution is stationary and no currents exist. 

Next, let us turn to non- equilibrium steady states. The simplest example is a system in contact with two energy 
reservoirs, e.g., two thermal baths at different temperatures: T\ and T^. For physically realizable cases, we can think 
of a time frame during which energy flows through our system (from the hotter to the cooler bath) steadily, so that 
the average energy within the system is time- independent. Typically, to model the coupling of our system to such 
reservoirs, there is no need to respect detailed balance when specifying a set of rates. Of course, the rates are not 
completely free of contraints. Without going into details, let us simply restrict our considerations to those rates which 
eventually take the system to a unique time- independent distribution: P*(C). However, without detailed balance, we 
should expect, generically, P*(C) ^ P* q (C). Now, a time-independent distribution does not imply that all currents are 
zero. Instead, we can expect (time independent) current loops to be present. Using the analog above, this situation 
corresponds to magneto statics, in which the charge density is stationary but steady current loops prevail. In this 
sense, we prefer to use the term "steady states" for describing such (magnetostatic-like) systems, while reserving 
"stationary states" for those in equilibrium (or electrostatics) . The presence of current loops also highlights another 
crucial difference between equilibrium and non-equilibrium steady states, namely, time reversal invariance. Since 
currents change sign under time reversal, only states with zero currents are invariant. In this light, the concepts of 
equilibrium, detailed balance, and time reversal are intimately connected. 

Once P*{C) is known, we can compute the "internal energy" of the system 

u = {H[c}) = J2n[c]p*[c}, (ii) 

c 

since we have a well defined microscopic Hamiltonian. For equilibrium states with P^ q {C) = e~P n ^ / X)c e~ l3n ^ c \ 
it is a simple step to check that, without knowing the details of TC[C], the specific heat: dU/dT oc —dlf/d(3 oc 

l^{TL[C] — U) 2 ^ is never negative. Notice that dU/dT > is true for systems of any size, regardless of whether a 

proper thermodynamic limit exists or not. For non-equilibrium steady states in, e.g., systems coupled to two baths, 
we naturally have U{T\,T2) and may extend the usual definition of specific heat. Of course, in this case, there are 
two such response functions: 

as we can separately "dial up" the temperature of either bath while keeping the other fixed. While we cannot expect 
both to be negative, we will show that one of them can be negative in very simple systems. 



III. SIMPLE SYSTEMS WITH NEGATIVE SPECIFIC HEAT 



A. A minimal system with three states 

We begin by considering an abstract system with only three, non-degenerate, energy levels. For reasons shown 
at the end of the next subsection, it is impossible to construct a system with only two microstates which displays 
negative specific heat. In this sense, we believe that the three state system is "minimal." 

Each configuration (or microstate) is associated with a unique energy. Using the subscript a = 0, 1, 2, let us denote 
the energies by 
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E a — 0, ex, £2 



(13) 



with £2 > £1 > 0. Our goal is, given a set of rates R(a 



to find U, the average energy in the steady state. First, 



we must compute P*, the probability for finding the system in level a when steady state has been reached. Then 

U = E a P* = Sl P? + £ 2 P 2 * . (14) 
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Had we been interested in the equilibrium distribution, corresponding to our system being in contact with a single 
thermal bath, then we could use, e.g., the usual Metropolis rates (Eq. ||). Here, they are explicitly 
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where <5 = £2 — £i- The reader may verify that, with these rates in Eq. (|J), the Boltzmann form, Pj* 2 = Pq exp (— /3e Q ), 
is indeed a solution to L \P*) = 0. Note further that, even if we forbid the <-> 1 transition (by setting R(0 — > 1) = 
i?(l — > 0) = 0), the system will still equilibrate to the same set of P*'s. 

Next, let us modify the dynamics by coupling the 0^2 and 1 <-> 2 transitions to different baths, while forbidding 
the «-» 1 transition entirely. To avoid confusion with subscripts, we will label the two bath temperatures by T x and 
T y . Defining 



x = exp [-£2/ / k B T x ] and y = exp [-8/k B T y ] 
the master equation takes the simple form 



d t Po 
d t Pi 



-xP 

-yPi 



P2 
P2 



d t P 2 = xP + yP x 

The steady state distribution is trivial to find: 

P£=y/Z; Pt=x/Z; P 2 * 
where Z = x + y + xy. Thus, the average energy is 

_ xsi + xye 2 
x + y + xy 

so that the specific heats (associated with the x- and j/-baths) are: 



2Po 



xy/Z 



C x 



dU 
df T , 



{ei + y£ 2 } 



xye 2 



and 



dU 

C u=ffT = i X£2 " £ i (! + x )} 



k B Tl [x + y + xy] 
xyS 



k B T 2 [x + y + xy\ 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



While the first of these never goes negative, the second clearly becomes negative for a range of x. Working out the 
details, we find C y < (for all T y \), provided T x drops below the critical value 



£2 

k B 



In 



£2 



£1 



(22) 



Note that, unless baths with negative temperatures are invoked, £2 > 2e\ is needed. 

While it is easy to analyz Eq. (21) mathematically, it may be helpful to provide an explicit example. We choose 
£2 = 4£i and plot, in Fig. 1, C y /k B against T y (in units £\/k B ) for various TVs. Note that, in this "minimal" system, 
C y does not change sign as long as T x is fixed. 
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FIG. 1. Specific heat C y /k B vs. T y for various T x 's. Starting from the lowest curve, 
fesT^/ei = 1, 2, 3, 4 respectively. The heavy line at the top is associated with T x = oo. 

B. An intuitive picture 

To help the reader appreciate how such an unusual response arises, we provide an intuitive picture for this phe- 
nomenon. Since the effects of a thermal bath on the transition between two levels are written in terms of rates, we 
may think of the role of a bath as that of (a pair of) "pumps" on the "population" in these levels. If we had only two 
levels, there could be only one pair of rates (pumps), regardless of how many "baths" are coupled to them. In the 
steady state, only one ratio of populations (say, the occupancy of the higher level to that of the lower) is relevant. 
In this situation, the only sensible way to define "thermodynamics" of the system is to have a single bath with an 
effective temperature. Then, the usual intuitive picture easily emerges: if this (effective) temperature is increased, 
the relative occupancy of the higher level will be higher, leading to a positive specific heat. The steady state system 
might as well be described as "in equilibrium." Therefore, a minimal system which can display negative specific heat 
must involve more than two (distinct) levels. 

With three levels, there are three (sets of) rates, so that our system may be coupled to three different "thermal 
baths." By forbidding the transition between the lower pair of levels (as in the simple case above), it is already possible 
to generate "unconventional" population changes. This restriction also reduces our consideration to only two baths. 
Now, let us assume that our system has settled in a steady state where the populations are given by (Po, Pi, P2), using 
the notation above. If we increase the temperature of the bath coupled to 1 <-> 2, we will certainly deplete Pi in favor 
of P2. However, the 0^2 transition is coupled to a different bath, with a fixed temperature, so that any increase 
in Pi will be driven to increases in Pq. If this bath is very cold, then the downward pump is very strong (relative to 
the "up-pump"), ensuring increases in P2 be accompanied by much larger increases in Pq. Thus, there is an effective 
"downward" shift of population (from Pi to Po), even though the thermostat coupled to 1 <-> 2 has been turned up! 
By contrast, if both transitions are coupled to the same bath, then there would be also an "upward" shift (from Po 
to P2), leading to the expected increase in U. From these considerations, we see that forbidding the 1 <-> transition 
introduces an "barrier" (or "obstruction" ) , leading to a "local minimum" . The presence of such "obstructions" forms 
the key ingredient to the display of negative specific heats. In the next section, we provide another illustration of 
such a link, in the context of currents and drives. Of course, "the proof is in the pudding." The heuristic argument, 
hopefully, helps the reader to digest the mathematical "pudding" above. 

Before turning to the explicit example, we note that this simple three state system can be generalized to include 
a third bath, coupling to the <-> 1 transition. To be specific, let the temperature be T w and, in addition, let the 
"efficiency" of this coupling be 77, so that the extra rates are 

R{0^l)=r)w; R(l 0) = r\ (23) 

with 

w = exp [-ei/feTj . (24) 

Of course, setting 77 to zero will reduce this system to the simple one above. Computing the steady state distribution 
and the average internal energy is a straightforward exercise. Here, we only quote the relevant result. For sufficiently 
small 77, there is a critical line in the (T x -T w ) plane, given by 
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on which the specific heat C y vanishes. For temperatures (or x,w) below this line, C y is negative. 



(25) 



C. An explicit example: the Driven Lattice Gas 

A "realization" of the above three-level system is the "two-temperature Ising model," || which consists of a lattice 
gas coupled to two baths. The only extra complication comes from degeneracies. Otherwise, all analysis is identical 
to the minimal case. 

Our system consists of a 2 x 3 lattice, half filled with particles. Imposing the boundary conditions: 

brick wall in x ; periodic in y 
the 20 possible configurations fall into three distinct classes, labeled as follows: 
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For energetics, we use the usual ferromagnetic Ising Hamiltonian: —J is associated with each nearest neighbor 
particle pair. For convenience, shift all energies by — 3J. As a result, the labels of the configurations are also the 
energy levels (in units of J). The degeneracies are, respectively, (2, 12,6) due to translation and parity (in x and y). 
To be pedantic, let us use a second label to distinguish them, so that all 20 configurations are, explicitly: 
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Next, let us specify the dynamics: particles may hop to nearest neighbor empty sites (Kawasaki exchange |4[|) 
according to the usual Metropolis rates (Eq. ^|). The two temperature model is defined so that exchanges along (x,y) 
axes are coupled to baths with temperatures (T x ,T y ). As examples, the rate for 0,1— >2,1 is exp(— 2 J/kT x ) but the 
rate for 1,1 — >2,1 is exp(—J/kT y ). Obviously, our system is reduced to the trivial equilibrium case when T x = T y . 
With these rates, the master equations can be easily written. Using subscripts to denote configurations, we give only 
three of the 20 master equations: 

dtPo.x = P2.1 + P2.2 + P 2 ,3 - 3xPo,i 

dtPx.x = P2.X + P 2 ,3 + Pi,2 + Pi, 7 + Pxa ' (2y + 3) Px.x (26) 
d t P2,x = xPq,i + V (Pi,i + Pi,2 + Px, 3 + Pi.e) + P2,5 + P2.6 - 7P 2 ,x 

where 
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x = e~xjp(—2J/kT x ) , y = exp(— J/kT y 



(27) 



To find the steady state P*'s, we note that all degenerate configurations will have same probability (since they are 
related to each other by symmetries). So, we just write equations without the second subscript: 



= 3P 2 * 

= 2p: 



3xP * 



= xP * + AyP? - 5P 2 * 
Apart from degeneracies, the solutions are same as in the minimal example above: 

p o = y/ z ; p \ = xZ ; p 2 = x v z . 

with 



The average energy is given by 



so that 



Z = 2y + 12x + 6xy . 



U _\2x + 2x6xy _ ^ l + y 

J Z y + 6x + 3xy ' 



a 



dU 



y ~ dT u 



oc [3x — 1] 



(with a positive definite proportionality factor). Similar to the minimal case, we have 



C y < 



(28) 



(29) 



(30) 



(31) 



(32) 



(33) 



for all T y , provided T x is low enough. Note that the degeneracies of these three levels affect the critical T xc . A simple 
exercise including g a (the degeneracy of level a) in the computations above leads to the general result 
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90 



(34) 



Interestingly, this critical value is independent of the degeneracy of the middle level. From the intuitive picture given 
above, this behavior is, arguably, reasonable. After all, the transition from positive to negative specific heat depends 
how the depopulation of the middle level is redistributed between the upper and lower levels. At the critical point, 
this redistribution leads to an average energy (associated with the outer levels alone) which is identical to Ex, leading 
to no change in the total energy. Now, each state in the middle level would "suffer the same fate," so that it is 
irrelevant how many such states there are. Thus, we have g\ independence. 



D. Energy flux 

In general, a system coupled to two thermal baths will lead to a transfer of energy from the hotter bath to the 
colder one. On the other hand, by definition of steady state, the energy stored in our system is a constant. Thus, we 
may expect a constant energy flux through our system. Most remarkably, in our minimal (77 = 0) case, there is no 
energy flux through the system! In this sense, our system and both baths are "in equilibrium," and detailed balance 
is satisfied! Indeed, the general condition for detailed balance [pj is less restrictive than Eq. (|7|), so that it is possible 
to have a system coupled to two baths and retain detailed balance. The ramifications of such "effectively equilibrium" 
systems remain to be explored. 

Meanwhile, as T y is increased (with fixed T x < T xc ), the energy of the system decreases. To be more precise, 
suppose we take a system in steady state and suddenly replace the y-bath by another with a higher temperature. 
Then there must be a net flow of energy out o/our system somehow, since the new steady state must have lower U. 
This paradox cannot be solved by studying the properties of steady states alone. The resolution lies in an analysis of 
the full dynamics. Deferring all details to another article, we provide only the result here: less energy is transferred 
from the y-bath to our system than from the system to the x-bath. In the general case (rj > 0), there is a net energy 
flux through our system, so that similar surprising behavior appears less paradoxical. 
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IV. A MODEL WITH NEGATIVE CONDUCTIVITY 



The counter-intuitive phenomenon of negative response is not restricted to the energy-temperature variables. In 
this section, we show a similar behavior in the current-drive variables, i.e., the possibility of "negative conductivity". 



A. Set up and the unobstructed case 

Consider a single particle hopping in two "lanes" of L sites, periodic on the long side and "brick walls" bounding 
the lanes. For convenience, assume L is even, so that the sites along the lanes can be labeled by the integer i = 
—L/2, . . . , L/2 — 1. Periodicity allows us to use i = 0, 1, . . . , L — 1 interchangeably. 

Let us denote the probabilities of finding the particle on the two lanes by 



P r and Qi 



Symbolically, we represent this system by: 



0-2 


0-1 


Qo 


Qi 
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P-2 


P-l 


Po 


Pi 


P2 



(35) 



(36) 



Let there be no obstructions and suppose the particle is "driven" by a uniform external "electric" field E, so that 
the jump rates are proportional to 



1 for crossing lanes only; 
x for jumping "upstream" (i - 
1/x for jumps "downstream" (i 



i — 1) to either lane; 
> i + 1) to either lane. 



(37) 



The normalization factor is just 



N = 



1 + 2x + 2/x 



We may regard a; as a Boltzmann-like rate: 



x = e 



-0E 



(38) 



(39) 



in which case N = 1/(1 + 4 cosh (3E). Note that the behavior of this "charged" particle is not the same as one in 
free space. Instead, it represents an overdamped situation, where inertia can be neglected compared to (thermal) 
damping. As a result, the particle settles down to the drift velocity instantaneously. Furthermore, with the rates and 
normalization chosen, the limit of E — > 00 does not mean the particle will move with infinity velocity. From the rates, 
it should be clear that, in this limit, all that happens is that the particle moves to the next site (either lane) with 
unit probability. Thus, the velocity saturates at unity for "infinite drive". 

Without obstructions, the system is translationally invariant and the steady state is trivially given by 



P t = Qi = 1/2L 



(40) 



since the particle can be located at any of the 2L locations. The uniform current density, J/ ree , can be found by 
considering the probability of all jumps between i and i + Taking into account the contributions from both lanes, 
we have -N {P t + Qi) - 2xN (P l+1 + Q l+1 ). Therefore we obtain 



Jfr 



2 sinh j3E 



(1 + 4cosh/3S)L 



Needless to say, the conductivity is positive 



dJ_ 
~dE 



> 



(41) 



(42) 



for all finite f3E. 

Another perspective is to consider the average velocity of the particle. Since we have only a single particle, we have 
J = pv. Given uniform density (1/2L), we obtain v = 4 sinh (3E/{1 + 4cosh/3i?), which is a monotonically increasing 
function of the drive E. Note that this expression is expected for both small and large E. In the former case, four 
out of the five possible jumps contribute to the velocity, so that the conductivity is simply (4/5)/?. 
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B. Single obstruction case 



Next, let us introduce barriers at i = 0, but only for one lane. Specifically, "impenetrable walls" are placed between 
the site associated with Pq and three of its neighbors (those associated with Qq, Pi, and Qi). Symbolically, we draw 
thick lines around Pq (Fig. 2), and, using the analogy of water flowing under gravity, we will refer to the obstruction 
as a "cup." Clearly, the cup causes a blockage, which, for large E, will reduce the current seriously. However, as E 
is lowered, the particle's chances of backward jumps are higher and leads to a higher current, giving us a negative 
conductivity. 
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FIG. 2. Schematic display of the probability for finding the particle in each box. Thick 
lines borders are associated with impenetrable walls. 

To see that this picture is indeed correct, we start with a careful consideration of all the master equations. Without 
the cup, these equations read, for all i, 



d t P t 



1 



{Pi-i + Qi-i) + x {P l+ i + Qi+i) + ' 



i 



2 

x 



2x P 



(43) 



d t Q l = ^ (P-i + Qi-i) + x {Pi+i + Qi+i) + Pi - ( 1 + | + 2x ) Q-, 



With the cup in place, the only jumps affected are those involving <-> 1. As a result, the only equations affected are: 



d t P = - (P-i + Q-i) - 2xP 
x 

d t Q = - (P-i + Q-i) +x(P 1 + Q 1 )-2[- + , | Q {] 

x \x 

dtPl = -Qq + x (P 2 + Q 2 ) + Qi- (l + -+x) P x 

x \ x 

d t Qi - -Qo + x (P 2 + Q 2 ) + P ± - ( 1 + - + x ) Q x 

X \ X 



(44) 



A factor of has been absorbed into the time scale, an irrelevant concern for the time-independent state. Setting 
the left hand side of these equations to zero, we seek the unique solution corresponding to the steady state. Note 
that, if we set the drive to zero (x — ► 1), then the obstruction becomes irrelevant in steady state and we retrieve 
the expected result: the trivial flat distribution. For the E case, we defer all details of the calculation to the 
appendix, providing a few brief steps and the result here. 
For convenience, define the sums and differences: 



S = P + Q 
D = P-Q 



(45) 
(46) 



With no obstructions, it is intuitively clear (and mathematically easy to show) that the two lanes have equal proba- 
bilities in the steady state. However, the cup enhances Pq, so that 



A = for 
l-x 2 



Do 



1 + 2.T 2 



So > 



(47) 



Invoking (probability) current conservation, we set 



K = -Si-!-2xSi i^0,l 
x 



(48) 



9 



which is proportional to the current between i — 1 and i. But this K must also hold for all links. In particular, 



K = -S-i - 2xS 

x 

is 2 r\ a So— D 
K = -Qo - xbi = 

x x 



(49) 



xS\ 



Using Eq. 



and 



?|), we find a set of equations for the S"s, leading to 

(1 + 2x 2 ) (2 



S = 



K = 



3Lx 2 + 2 + 2x 2 ~ Lx 2L - 2Lx 2 + 2L - 2x 2L - 2x 2 + 2L 

2x(l-x 2 ) (3 - x 2 ^- 1 ) - 2x 2L ) 

3Lx 2 - Lx 2L - 2Lx 2+2L - 2x 2L + 2- 2x 2+2L + 2x 2 



(50) 



(51) 



It is instructive to check that these expressions indeed reduce to 1/L and 0, respectively, in the limit E — > (x — > 1). 
Another interesting limit is f3E diverges faster than \nL (Lx — > 0). In this case, the probability of the particle jumping 
out of the cup is so low that even the entropy factor (the 2L possible locations for the particle) is not enough to 
overcome the barrier (x -t 0), so that "complete trapping" occurs. This expectation is borne out by 



and, from Eq. @, [D 



trapped 



[So] 

trapped 

1 also. Thus, we have 



[Po. 



trapped 



1 and [Qo] 



trapped 



o. 



(52) 



(53) 



Verifying that all other probabilities also vanishes, we see that the particle is completely trapped in the cup. To 
complete the study of this limit, we find 



K tr 



ed = 3x + 0(Lx 2 ,x 2 ) -> 



(54) 



Since J = NK, the current vanishes (as x 2 ) for large fields, rather than saturating, as in the free case. Thus, d,J/dE 
must be negative in some range of E. 

Returning to Eqs. (§Cj) and (fl|), we find less cumbersome forms if we consider the thermodynamic limit and drop 
terms of 0(x 2L , Lx 2L ) first. The results are: 



So 



(l + 2x 2 ) (2-x 2 ) 



3Lx 2 

Jobstruction NK 



2x 2 



6x 2 (l-x 2 ) 



(2 + 2a; 2 + 3Lx 2 ) (2 



2x 2 ) 



(55) 
(56) 



0.05 




FIG. 3. Current vs. i: 



-HE 



Note increasing E corresponds to decreasing x and the 



linear rise of the current for small E (near x = 1). For large enough drive (small x), the 
current drops with increasing E. 
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Of course, we can differentiate the current and find the critical field above which the conductivity dJ/dE is negative. 
However, it may be more instructive to provide an illustration of this expression. In Fig. 3, we plotted this current 
as a function of the driving field, for the case of L — 10. We should remind the reader that x — exp(—/3E), so that 
larger x corresponds to smaller E. Note that the curve is basically linear near the x = 1 end, confirming the typical 
constant, positive conductivity for small E. On the other hand, negative conductivity sets in for sufficiently large E 
(approximately fc^T in this case). 

V. CONCLUDING REMARKS 

In this paper, we have shown that "negative responses," in a system coupled to thermal baths, can be easily 
generated, provided we consider non-equilibrium steady states. Specifically, we gave examples of exceedingly simple 
systems which display negative specific heat. One is an abstract system with only three energy levels, but coupled to 
two thermal baths. The other is a standard Ising lattice gas (3 particles in 2 x 3 sites), also coupled to two baths. 
These models are "minimal" in the sense that systems with just two energy levels can support only a single (set of) 
transition rates, so that an single (effective) temperature can be defined. Then, we have essentially an equilibrium 
system, with non-negative specific heats. There is no doubt that such systems can be extended to more complex and 
macroscopic cases. In particular, using Monte Carlo simulation techniques with the two-temperature Ising model, we 
have observed negative specific heats in, e.g., a 60 x 60 lattice for T y > 1.5Tq with T x = 0.7Tq- Another example is 
a two-temperature Ising model with a population of mobile impurities that couple to a heat bath with a temperature 
different from that of the Ising spins. This model has been proposed to mimic the non-equilibrium behavior of a 
biological membrane in which diffusing protein particles receive energy from a heat bath with a temperature that is 
different from the membrane matrix in which diffuse. Under certain conditions this has shown to lead to a negative 
specific heat (Fig. 4 in |(j)). We also showed a second form of negative response: negative conductivity. The explicit 
example consists of particles hopping along two lanes of discrete sites, driven by an external field so that a non- 
vanishing DC current is established in the steady state. When an obstruction is introduced, the current is shown to 
decrease when the drive is increased (beyond some critical value). In both cases, the key ingredient (besides being 
in non-equilibrium steady states) is the presence of local minima or "obstructions." When the drive, be it another 
thermal bath or an external field, is too large, the system is caught for longer in the obstruction. By lowering the 
drive, the system can respond "more positively". Though the phenomenon has been observed for some time, e.g., 
in driven diffusive systems , it has been discovered in another context and popularized through the catchy phrase: 
"freezing by heating" ||. Since this kind of counter-intuitive response is not only reserved for negative specific heat, 
but also negative conductivity and beyond, we may label all of them by the catch-all phrase: "getting more by pushing 
less." 

Though we are certain that this kind of phenomena are present in physical systems, we have not conducted 
a systematic search for them, nor have we proposed designs for such devices. The goal of this paper is simply to 
encourage those who teach statistical mechanics to keep their students' minds open about the issue of negative specific 
heats. In particular, in non-equilibrium steady states, positive responses are not intimately related to stability. 

To close, we should mention that "negative specific heat" is not a novel phrase. Indeed, it occurs frequently in 
discussions of self-gravitating systems. Evidently, Eddington wrote about how a star or star-cluster would cool down 
if energy is added |9). Reading the literature, there is no doubt that this is a subtle situation of a complex system, 
with particles subjected to long range interactions (gravity). Furthermore, it is clear that such systems are inherently 
unstable, leading to "gravothermal catastrophe" and collapse. Clearly, the theoretical existence of such a kind of 
negative specific heat was not well received by physicists. In the opening sentence of a recent article |U^, Lynden-Bell 
wrote, "When I first used the concept of Negative Specific Heat ...the Statistical Mechanics community thought I 
was talking nonsense." By contrast, our examples are far less exotic and the phenomenon should be abundant in 
our neighborhood as well as amongst the distant stars. After all, every living organism can be considered as a non- 
equilibrium steady state, coupled to more than one reservoir of energy so that there is a constant through-flux. If we 
succeed in convincing other teachers to add a footnote in their course, that dU/dT > is ironclad only for system in 
equilibrium, then this paper would have served its purpose. 

VI. APPENDIX 

Here we provide a few more details for finding the steady state distribution in the two-lane hopping model. Referring 
to the schematic diagram ([56]), we seek solutions to the set of equations (E|,Q) with zero on the right hand sides. 
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Recall that we have periodic boundary conditions, so that hopping rates from sties i = L — 1 to i = — L are the same 
as all others (except for the obstruction). 

As mentioned in the text, it is convenient to define sums and differences: 

S = P + Q 

D = P-Q (57) 



so that we have 



= + 2xS i+1 - 2 ( - + x ) St (58) 



x \x 



-2[l + -+x)D i (59) 



x 



and 

= -S-i + xSi - 2xS --(S - Do) (60) 

X X 

= -xSi - 2xD a + -(S Q - D a ) (61) 

x 

0=-{S -D ) + 2xS 2 -(-+x]s 1 (62) 
X \x J 

= - ^2 + | + xj D 1 (63) 

from Eqs. (^3|,Q), respectively. The advantage of this decomposition is obvious now. It proves the intuitive notion 
that, apart from the obstructed sites, the two lanes should have equal probabilities in the steady state. So, 

Di = for i ^ (64) 

Do = (65) 

Now, we have a set of equations involving only S: 



= —S-i + xSi - f + ^ X xSq (66) 
x 1 + 2x z 



= So + 2a;^ 2 - f I + x ) .S', ( 07 ! 



x 



= -5i_i + 2.T5 i+ i -2 f i+x) S* 4 for i ^ 0, 1 (68) 
x \x ) 

To solve these, we exploit current conservation, since the physical content of these equations lies in the difference 
between the particle currents into and out-of site i must be zero. Thus, we first consider the "unaffected" sites, just 
to be careful. Set 2Si-\/x — 2xSi to a constant: 

K = -Si-i - 2xSi i ^ 0, 1 (69) 
x 

which will be proportional to the steady state current. (As a reminder, given the periodic condition, Sl/ 2 = 5 , -l/2, 
K = 2Sl/2-i ~ 2xS'_ i /2 is P ar t of this set of equations.) Note that this "trick" solves Eq. (68) automatically: 
= K — K. Next, this K serves as the steady state current, so that it must hold for all jumps i i + 1. In particular, 
the jumps across —1^0 leads to 

K = —S-i — 2xS (70) 
x 

which is the same as the unaffected case. Meanwhile, the jumps across «-> 1 provides a new equation: 

K = -Q - xSi (71) 

x 
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By definition, 2Qo = So — Do, so that, using Eq. (|47|), we obtain 

3x 2 



2Qo 



l + 2x 2 



Sq 



and 



K = 



3.7' 



l + 2x 2 



Sq — xS\ 



(72) 



(73) 



All these equations with K on the left can be recast as recursion relations for the Si's. Specifically, we will try to 
express everything in terms of Sq, so that we write 



Si = j-^So - -K 
1 + 2x z x 

Sj_i = -if + x 2 S l for i^l 
Starting with the i = recursion relation Eq. (|75|), we work all the way around back to Si: 

S_i = | if + x 2 S Q 

S_ 2 = | K + x 2 S^ = ^(l + x 2 )K + x^So 



(74) 
(75) 

(76) 



S_ri 



| (l + ^^ + x^-^K + x^So 
x (1 - x 2n ) 



2(1 -x 2 ) 



■K + x 2n S 



Si — S-L+l 



,2(L-1)) 



-K + x 2 ^So 



2(1 -x 2 ) 

Together with Eq. ((74J), the last of these allows us to express K in terms of So- 



(1 



r 2(L-l) 



2(l-x 2 ) 



1 



K = x 



l + 2x 2 



x 



2(L-1) 



So 



(77) 



Finally, normalization (Y]. Si — 1) will fix everything. 

It may be instructive to compare the above equations with those in the free case. The only difference is that we 
would have one more "unaffected" step. Instead of Eq. (f74|), we have [So]^ ree = ^Kf ree + x 2 [Si]j ree , so that 



1 _ x(l-x 2 ( L ~^) 



2{l-X 2 ) 



+ - 2(i - 1} [So] fre 



i.e.. 



[3d. 



free 2(l-x 2 ) 



K 



free 



(78) 



(79) 



Note first that this expression is independent of L\ This result should be expected, since we know that the average 
velocity must be L-independcnt for the free case. But the velocity is just the current divided by the particle density 
(i.e., probability for finding the particle at, say, site 0). To continue the check, insert this expression into the recursion 
relations and find that all the S's are the same. So, normalization gives [Si]j ree — l/L and Kf ree = (4sinh/3_E) /L. 
Armed with these considerations, we return to the case with barrier and define 



M = So 



2{l-x 2 ) 



K 



(80) 
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which is a measure of the effect of the cup. Eliminating K in its favor, Eq. ( f77[ ) is replaced by a neater formula: 

(l-x 2L ) 



(1-x 2 ) 



1 



2 (1 + x 2 ) 



(81) 



Finally, the sum over the S's is 



L-l 



L-l 



71 = 1 



71=1 



2n \ 



X (1 — X 

2(1 -x 2 ) 



■K + x 2n S 



L-l 



2{i-x 2 y 



Eliminating K and setting this to unity, we have 

1 = LS 



2L 



1 — X 
l-x 2 



L 



M 



(82) 



(83) 



Eqs. ( pi| , p3| ) can now be used to obtain Eq. (50). 
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